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Abstract 



S 

-)— > 

a 

O I The Xy model (s — 1/2) on the one-dimensional alternating superlattice (closed chain) is 

solved exactly by using a generalized Jordan- Wigner transformation and the Green function 

method. Closed expressions are obtained for the excitation spectrum, the internal energy, the 

^ \ specific heat, the average magnetization per site, the static susceptibility, x^^i^nd the two-spin 

lO ' correlation function in the field direction at arbitrary temperature. At T = 0, it is shown that 

^1 ' the system presents multiple second order phase transitions induced by the transverse field, 

^^ , which are associated to the zero energy mode with wave number equal to or vr. It is also 

shown that the average magnetization as a function of the field presents, alternately, regions 

of plateaux (disordered phases) and regions of variable magnetization (ordered phases). The 

static correlation function presents an oscillating behaviour in the ordered phase and its period 

2 , goes to infinity at the critical point. 



2 ■ 1 Introduction 

O 

J-^ [ The experimental development of magnetic superlattices, by using molecular beam epitaxy tech- 

nique [1-3], has increased the interest in the study of these systems. Although they are three- 
^^ ■ dimensional systems, there is a predominance of the one-dimensional behaviour in their proper- 

j^ ■ ties, and this is the main reason for studying one-dimensional superlattices. Therefore, interest 

has been considerably increased in the study of spin systems on these lattices. 

Among the spin systems the XY-model (s=l/2), introduced by Lieb et al. [4], occupies a 
special place, since it can be solved exactly for the homogenous lattice. Although almost all 
static and dynamical properties are known for the model on the homogeneous lattice (see [5] and 
references therein), the known results for non-homogeneous periodic one-dimensional systems are 
restricted to the alternating chain [6-8] and to the excitation spectrum of the general alternating 
superlattice [9], and its critical behaviour, which has been obtained by using the position space 
renormalization group approach [10]. 



In this paper we consider the isotropic XY-model in a transverse field on the one-dimensional 
alternating superlattice (closed chain). We solve the model by introducing a generalized Jordan- 
Wigner transformation [9] and by using the Green function equation of motion tecnhique. 

In section 2 we determine the relevant Green functions and present a detailed discusssion of 
the excitation spectrum. In section 3 we obtain the internal energy and the specific heat. The 
induced magnetization is studied in section 4, and in section 5 we calculate the two-spin correlation 
function. Finally, in section 6, we summarize the results and present the main conclusions. 



2 The Excitation Spectrum 

The superlattice that we are going to consider consists of cells composed of two subcells A and B 
with UA and ub sites, respectively. The l-th unit cell is shown in Fig. 1. The distance s between 
two consecutive sites is taken as unity. 

If we assume periodic boundary conditions for a chain with TV" cells, the Hamiltonian of the 
XY model [4] can be written in the form 

l^\ \ m—l ni—l 



UA 

E 

m— 1 



'^hASl^yyi + 2Z '^^bSJ^^^ 



(1) 



where I identifies the cell, S^ = S^ zt iS'^ , J is the exchange parameter between spins at the 
interfaces, Ja{Jb) the exchange parameter between spins within the subcell A{B), and hAihs) 
is the transverse field within the subcell A{B). The spin operators can be expressed in terms of 
fermion operators using the generalized Jordan- Wigner transformation [9] and, by introducing this 
transformation, the Hamiltonian can be written in the form 

N (riA-l 



/^l l^ m—l m—l 

+ J [A,nA^l^ + ^lnBa' + l,l) + h.C. + 



E 2^-4 [al^ai,^ - 1) + £ 2/15 {hL^m - 1 



m—l 

where a's and 6's are fermion operators, and $, given by 

N UA 
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(3) 



is a boundary term which will be neglected. As it has been shown [11], this boundary term, in the 



thermodynamic limit, does not affect the excitation spectrum, the static properties of the system 
nor the dynamic correlation function in the field direction. Introducing the Fourier transforms [9], 



agfci = \ Jfl -—Y^^exp{iQdl)sm{mki)ai^m, (4) 

V l,m 



2 



^Qk2 



N{nB + l) 



— ^cxp(i(3(i/)sin(TOfc2)6i,m, (5) 



where ki = hitt/ (ua + 1) , ni = 1, 2, . . . , ha, k2 == n27r/ (ub + 1) , n2 = 1, 2, . . . , ub, Q == 2Tm/N, 
n — 1,2, ... ,N, d — HA + riB is the size of the cell, the Hamiltonian can be written as 

Q ^ ^ 

where Hq is given by 

fci fc2 

+ sin(fc2nB)sinA:i exp(i(5d)aj^j,^5(5fe2 + /i.e. j- , (7) 



j,A{B) _ 



with E^'^J^ = -Ja(b) cosfci(2) - hA(^B)- 

As in the study of the excitation spectrum [9] we will solve the model by using the Green 
function method [12]. Adopting the notation ((3fJi(i);3fJ2(0)))^ for the retarded anticommutator 
function, where Sfti and 5ft2 are arbitrary operators, and introducing the time Fourier transform 
defined as 

{(JRi(tyJMO)))r = ^ f ((3?i;3?2))cxp(-iL^t)du;, (8) 

2t^ J-oo 

we can write the equation of motion for the Green function ^ CLQki ; flo't' ^ in the form [12] 

- / ) [smA:2Sm(nAfci) + 

V(nA + l)(nB + 1) t^ 

+ (zQrf) sinfcisin(nBfc2)] (^{^6Qfe2;aQ,fe/ )^ , (9) 

where we have assumed fi = 1. 



Likewise we can write for <C &Q*:2; '^o'fei ^ ^^^ result 

((^Qfc2;aQ'fci)) = ,, '^^-.J' ==^[sinfc2siii(n^fc;) + 

+ exp{-iQd)smk[sm{nBk2)]((aQk'^;al^,f^)\. (10) 

Eqs. (11) and (12) constitute a closed set which can be easily solved by introducing the 
operators Aq „ and BQn defined by 



Aq^n ^ J — -^sin(nfci)aQfci, (H) 



ki 



BQ,n ^ J ^ ^ sin (nfcz) bQfc, (12) 



UB . - , 

fe2 



Therefore by eliminating the function ^ ^Qfe2; ^o'fc' ^ from eqs. (11) and (12) and introducing 
the operator ^q,„ we can find the result 

+ [fEl i^) Itn. i^) + cxp {iQd) f,% (co) fti [u^)] ((Aq,„,; 4,_„)) + 
+ [/m {^) fti (^) + exp {-^Qd) iLs (^) /™,n. (^)] ((^Q,i; 4'.")) 



where 



(13) 



.A / N _ J .^ sm(TOfci)sm(nfci) 

y/iuA + 1) (ns + 1) ^ ^ - £^fci 

.5 / X _ J_ >^ sin(TOfc2)sin(nfc2) 

^{nA + 1) (ns + 1) ^ w - ^fc. 



From eq. (15) we can find immediately that ^ AQ^rn',AQ, „ ^ is given by 



.t \\ _ 2 /ns + 1 



\ —Th,Q' {fm,n ^ + p f N X 

X { [fi,i H fta. {^) + exp {iQd) A^„^ {lo) fti (^)] X 

X [ft H {ft A i^) fli (^) + ft (^) Its (^) c^P {-^Qd)) + 

+ /n^^n. (^) (1 - ft (^) /m (^) - /l^n. (^) /m.. (^) ^^P (-^0^))] + 



[/:^- 



M (^) /™,i (^) + exp (-iQd) ft, (u;) f^n. (u;)] x 



X [ft. (^) (ft A (^) /m (^) + ft (^) /l^'nB (^) C'^P i'Q'^)) + 

+ ft (^) (1 - ft (^) /m (^) - ft. (^) .C« (^) exp Wd))] } } 



(16) 



The details of the calculation can be found in ref. [13]. 

By a similar procedure we can write the set of equations 



^ V" [sin fc2 sin (n^fci) + 

^{UA + l){nB + 1) V 

+ exp(-«(5rf)sinfcisin(nBfc2)] //aQfei;6^,j,, \\ , (17) 



7 f TT'A \~^ 

'^Qki-^^Q't')) = n "^^ru' ==^[sin/c2sin(n^/ci) + 

+ exp(i(5(i)sin/cisin(nBfc2)] U 6Qfc2;&Q/fcW h (18) 

which can be solved by eliminating the function ^ agfei; ^g/j.' ^ in the previous equation and by 



introducing the operator i?Q^„. This procedure leads to the set of equations 



BQ,ni]BQ, ^ 



2 luA+l B / ^r 

= -f\ —rJm.n{^)^Q,Q' + 

J \ UB + 1 

+ [ft (^) fins H + exp (iQd) ft. {u) fl, {lo)] ((Bq,„«; B^,^„ 
+ [ft (^) /m,i (^) + exp i-^Qd) ft. (u;) C„^ (u;)] ((Sq,i; ^S'," 



(19) 



and, as we can see, it can be obtained from eq. (15 ), provided we make the substitution A^-B 
and B^-A. Therefore ^ Bq^rn^BQ, ^^ ^ is obtained from eq. (18) by introducing the previous 
transformation. 

The excitation spectrum is given by the poles of ^ Ag,™; Aq, „ :^, or ^ Bq^rn', Bq, „ :», and 
corresponds to the solution of the equation Rq (lu) — 0, which is explicitly given by [9] 

1 - ^ft. (^) ftn (^) COS (Qd) - 2ft (io) ft H + 

+ [{ft (-'))' - (ft. i^)f] [(.ft i^)f ifts i^)T] = 0- (20) 

As mentioned in our previous paper [9], it should be noted that the values uj = E^ and 
u! — E^^ are not poles of <C Aq^m] A\^i „ ^, since in this limit the function is finite. Therefore, as 
expected, the spectrum of each subcell does not coincide with spectrum of the superlattice, and 
contains ua + ub branches [9] . 

The general solution of this equation is determined numerically, although analytical solutions 
can be found for some special cases. For instance for ua — nB ^ 2,hA = hB = h, we can express 
explicitly the solution in the form 

LUq^-h±^^c±^g{Qd), (21) 

where 

c - Y + \iJA + Jl), (22) 

giQd) = ^J^JaJb cos {Qd) + ^{JI + Jl) + ^{JI-Jl)\ (23) 

which for Ja = Jb reproduces the known results [6,7]. 

In the homogeneous medium limit, the spectrum obtained from the eq. (22) presents ua +nB 
branches, which correspond to a {ua + nB — l)-folding of the spectrum. As an example of this 
limit we present in Fig. 2 the excitation spectrum for ua = 2, hb — 3, Ja = Jb — J = ^ and 
hA^hB ^ 2. 

For zero field and ua + n-s odd, there is a zero energy mode with wave number different from 
zero, as shown in Fig. 3 for ua = 2, ub = 3, Ja = 2, Jb = 3 and J = 1. On the other hand. 



for ua + ns even this mode is not present, and for each wave number Q there are symmetrical 
solutions +ujQ and —ujq, as can be seen in Fig. 4 for ha = tlb = A, Ja = 0.3, Js = 3 and J = 1. 

As it has already been shown [14], the effect of a homogeneous field, Ha = hs ^ h,is to shift 
the zero field spectrum. This can also be shown directly from eq. (22) and, consequently, the 
existence of a mode of zero energy will depend on the strength of the field. 

As we can also see in Fig. 4, the extreme bands of the spectrum are very narrow, and this is 
related to the difference between the exchange parameters of subcells A and B. As the difference 
between the parameters of the media A and B decreases, the dispersion increases, and the gaps tend 
to zero. This is the expected behaviour, since we have a maximum dispersion in the homogeneous 
limit, as shown in Fig. 2. In all cases presented the number of branches is equal to the number of 
sites per cell. 

Again, we note that the spectrum can also be calculated exactly by using the position space 
renormalization group approach [10], and approximately by using a transfer matrix method [15]. 
Although the latter is an approximate method, we have shown that it reproduces numerically the 
exact result [5,8]. 

3 The Internal Energy and the Specific Heat 

The operator Hq defined in eq. (9) can be written in terms of operators Aq^n and Bg^n in the 
form 



Ho = 



UA 



riB 



J r 



jJl^k, 



jJl^k, 



^sin(nfci)A|5„ 

n=l 

J2 sin (nfca) B^_„ 



UA 



y^sin(nfci)^Q^„ 

n=l 
riB 

^sm{nk2)BQ^n 



.n=l 



Aq.ua^Q-^ + e^P (^Q^) ^qaBq-ub + h.^ 



(24) 



From this equation we can see that the internal energy, {H) —^n {Hq), can be calculated from 
the Green function Gq (w) given by 



Gq (lu) = Gig (uj) + G2q (lo) + G3q (lo) , 



(25) 



where 



G\q (c.) 



fci 



y^ ^sin(mfci)siii(nfci) ( Mq^,„; A|j^ 



G2,M - ;^^TT^^' 



^ ^ sin (mfcs) sin (nfca) {l^BQ,m;Bl^^^ 

rn—l n—1 



G3qH = -77[((SQi;^kn.))+CXp(iOd)((BQn,;A^, 



(26) 



(27) 



^QnA ; Sga)) + cxp (-^Qd) ((^Qi; ^^ 



(28) 



The Green functions <Bqi;AJj„^ >,<Bq„^;A^_i »,<ylQ„^;B^^ > and< Agi; B^ „^ > 
can, with the aid of eqs. (11-14) and (19,20), be written in terms of the functions ^ j4q ,„; aL, ^^ ^ 
and <C Bq^m', Bq,.^ >, and are given by 



^«,i;B|3., 



^[ftj.,{{B,,X...}} 



+ exp(^Qd)/i4l(c.)((i?Q,„,;i?1J,.„ 



(29) 



^Q,nA;^Q',l 



na + 1 r 
riA + l 



/mH((sq,i;sSm 



CXp (iQd) /i^„^ (w) {{BQ^nn;Bl^,,, 



(30) 



Bq.ubj ^Q',1 



riA + l 
nB + 1 



/f„«H((^Q,".;^^M 



+ exp(-zQd)A^i(L^)((ylQ,i;A|3, 1 



(31) 



■Sqi;^q',«^ 



;^[/f.M((A,„.;4^„.)> 



+ exp(-*Qrf)/5,^(a.)((AQ,i;4, 



(32) 



Then by using eqs. (18), and (31-34) in eqs. (28-29) we obtain the Green functions, GIq (oj) . 



G2q (lu) and G3q (uj) , which are given by [13] 



giqH = J2 



K 



Jc 



E' 



K 



X {sin^ fci [fl, {u) ft, H ((/fi {Lo)f [Its i^))' 
+ sin fci sin (n^fci) [li,nB i^) ^^^ (Q'^) + 

+ /i^„.H((/fiM)'-(.CH)')]}, 



(33) 



(/fiM)'-(/f„.H)']}, 



(34) 



and G2q (uj) can be obtained from GIq (w) by introducing the transformations ki ^ k2, A ^ B 
and B ^ A. From these results we can show that the expression for Gq {uj) has the form 



Then, the internal energy, (H) ~ J^o {^q) : can be obtained from the expression [7] 



(35) 



OO 



Im [Gq ju)) 
e/'^ + l 



du), 



(36) 



where, as usual, (3 — l/ksT. 
By using the Dirac identity 



r 



=F inS {x — xq) , 



X — xq ± ie \x — xo ^ 

we can show immediately that Im {Gq (uj)) can be written in the form 

Im {Gq {uj)) = ^ Fq^J {uj - UjQ^r) , 



(37) 



(38) 



where 



FQ,r — 



Rq (^Q,r) Gq {ujQ^r) 



(39) 



where r labels the branches of the spectrum and R'q {uj) is equal to (IRq {uj) /du 



The specific heat is obtained from eq. (38) and can be written in the form 

oo 

N dT Ntt ^ J (e/3- + 1)2 ^ ' 

The behaviour of the internal energy as a function of temperature does not present any re- 
markable difference as we change the superlattice parameters. A typical result is shown in Fig. 5 
where the internal energy is shown as a function of temperature T*{T* = ksT) for ua — nB — 10, 
Ja = '^■,Jb = '^, J = 1, hA = hs = 1.5. 

On the other hand, the behaviour of the specific heat as a function of temperature is very 
susceptible to these parameters. It can present a single peak or a double peak, as it can be 
seen in Figs. 6 and 7. This important feature, namely, the appearance of the double peak, is a 
consequence of the packing of the branches of the excitation spectrum, which is strongly dependent 
on the interaction parameters. As we increase the field, we move the spectrum downwards and 
this has the effect of suppressing one peak of the specific heat. This is shown in Fig. 8, where we 
have used the same lattice parameters of Fig. 7, and a larger field. 

In Figs. 9 and 10 we present the low temperature behaviour of the specific heat shown in Figs. 
7 and 8, respectively. For the results shown in Fig. 9, the excitation spectrum has a zero energy 
mode contrary to the case shown in Fig. 10, where it is not present, as can be verified by using the 
analytical solution presented in eq. (23). As expected, the derivative of the specific heat, dC/dT* , 
at T = 0, is equal to zero only when there is a zero mode energy in the spectrum. 



4 The Induced Magnetization and the Susceptibihty x^^ 

The local induced magnetization which corresponds to the average value of Sf'^ and Sf!^, is given 
by 

which can be determined from the Green functions ^ o.j.m',ajm ^ ^^^ ^ ^j,in',bj ^^^ 3>, respec- 
tively. In terms of the operators Aq^m and Bq„i these functions can be written in the form 

«^-;4™)) = 4 E expH(0-0')rf]((^Q,m;^lj',™)) (42) 



Q,Q' 



and 



(hm-^bl^)) = ^ E expH(Q- Q')rf] {{BQ,m;BlJ)) . (43) 



Q,Q' 
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Then introducing in eq. (44) the function <C Ag^^; AL ^ » , which is given in eq. (18), we can 
write ^ flj,™; a] ^ » in the form 



],m 



O-j.m'i '^j.m 






Rq{u 



[A^™ (^) /™,„. (^) (2,C (^) cos(Qd) + 2A^„^ (^) 



X ((A^, (c.))^ - (A^„, (-))')) + ((A^™ (-))' - {fin. (-))') 
X (A^i (-) + ft (-)) ((/fi (-))' - (/f„« (-))')] } . 



(44) 



The function <C &j,m; K ,„ :^ is obtained from <C Oj ,„; at ^ by introducing the transformations 
jA^^jB^jB^^jA^^^^^ "•B and TiB -^ riA- From this result we can write 



1 



Im ( {ai„,-a 



'3,"i' ^j,m 



dw. 



which, in the thermodynamic limit, allows us to obtain {S^^) from the equation 

F{q,ujq.r,h:) 



(45) 



^tm)- o + T^Y. 



2-K 



2 ■ 27r^7o (e/3-,.^ + 1) i^^M 



-dq. 



(46) 



where 



F{q.^q,r.h) = |y^£±I{A^„(c.)/^,„,(c.)[2A^„^(a;)cos(g) + 2A^„,(c^)x 

X ((A^i (-))' - (A^„. (-))')] + ((A^„ (c))' - (Al„. (-))' 

X (A^l (C^) + ft (-)) ((/m (-))' - {fir.. (-))')] } ' 



(47) 



and where q = Qd . The local induced magnetization in the subcell B, (S^^ ,is determined by 
following the same procedure and by using the function ^ ^j,m,\ bj„^ ^. 
The average induced magnetization in the cell, {S^^i), is defined as 



i^ceU — 



1 



N (ha + nB) 



N / 11/ 



Z^ Z^ \^l,m) + 2^ \^l,-m 



_l=l \m=l 



m—1 



(48) 



and from this expression, for the case Ha — hs — h, we obtain susceptibility x^^, which is given 

by 



X'' = ^ {Sl.i) 



(49) 
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Figure 11 shows the average magnetization in the cell and the susceptibility as a function of 
h at T = 0, for J = 1, Ja ~ '^, Jb — i, hA = hB = h and ua ~ «s ~ 2. As can be seen, the 
magnetization presents three plateaus which are limited by four critical fields, he, which correspond 
to the singular points of the susceptibility x^^ (x^^ -^ oo when h ^ he). These four singularities 
correspond to the modes of zero energy with wave number equal to or tt, which limit each energy 
band. For this case, the critical fields can be obtained exactly from eq. (23) and are given by 

h,e = ^^-0.691, 

/12c - ^^1|^- 1.096, 

V29 + 1 _ 
hzc = -. = 1.596, 

/i4c - t'^ - 1-809. (50) 

It should be noted that the plateaus, where the susceptibility goes to zero, correspond to the 
gaps in the spectrum. 

Figure 12 shows the average induced magnetization , at T = 0, for J = 1, J^ = 2, J^ = 3, 
hA = h-B — h and ua = 2, ne = 3, as a function of h. In this case, since the number of sites 
per cell is odd, there is a zero energy mode even for h = 0, and consequently there is no zero 
magnetization plateaux as in the case shown in Fig. 11. The critical behaviour is, naturally, also 
present when we have an inhomogeneous field. Ha ~ 1.25/ib = h, and this is shown in Fig. 13 for 
a larger unit cell, ua = 2, ns = 3. 

In Fig. 14 we present the magnetization for finite temperature (/?J — 50) for the superlattice 
considered in Fig. 13 with Ha = h-B = h. Although the central regions of each plateaux remain, 
as expected, the thermal fiuctuations suppress the quantum transitions. 

5 The Two-Spin Correlation Function in the Field Direction 

For sites in the subcell A, for example, the two-spin static correlation function in the field direction 
is defined by 



{^j,m^j+r,n ) ^ \^j.m^j-n^j+r,in^j+r,n) 



The average value (aj „jaj^„aL j,„aj_|_r,n) can be obtained from the expression 






duj. 



(52) 
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where the Fourier transform of the Green function <C o,j+r,n;aj^rn'^j,maj_^_^.„ '>, by using Wick's 
theorem, can be shown to be written in the form 

(a]+r,„aj,m) ((aj+r,n;a].„,)) ■ (53) 

Then, the correlation function ('S'^m'^j+r,n) can be obtained from the equation 



' —doj ■ 



e'^'^ + l 



'^j.m'^j-m} \'^j+r,n:0'j+r,n 



The Green function <C aj+r,n'.,aj^^ » can be obtained by using eqs. (6), (13) and (18), and is 
given by [13] 



ij+r.,„;4m)) = ]^Y^^7^IIc^PH90{/m,„('^) 



^{((/f„.(-))'-(/i"i(-))'lx 



X [ftn. {lo) /„^„, (C.) A^i {lo) - ft,n. {u) f,\ (c.) f,\^ (uj) + 
+ Itn H ft H ft (^) - fit (^) ft. (^) ft. i^)] + 

+ fts (^) [ft (^) /™,". (^) exp (^9) + /4.. (^) /i^m (^) exp (-zg)] + 

+ ft (^) [/m,n. (^) /«"!«. (^) + /l^™ (^) /i'n (^)] } } ■ (55) 

Therefore, introducing the previous result in eq. (56) and by using eq. (47) we can obtain 
numerically the direct static correlation function {Sf'^Sf^^.^.^) ~{Sf;'^) {Sf^^.^), and the results 
are shown in Figs. 15 and 16. 

In Fig. 15 we present the direct static correlation function as a function of r, at T = 0, for 

13 



two values of the field of the transverse field, h ~ 0.693 and h — 0.75, and for J = 1, J^ = 2, 
Jb = 'i, hj\ = fiB = h and ha — «b = 2, which are the same set of parameters of Fig. 12. The 
values of the critical fields, for this set of parameters, are given in eq. (52) and the first two values 
are h\c = 0.691 and /i2c — 1-096. As can be verified in Fig. 15, for these values of the field, the 
correlation function presents an oscillatory behaviour and its period increases ( period — > oo) as 
the field aproaches a critical value {h -^ he). 

In Fig. 16 it is also presented the direct static correlation function for the same set of lattice 
parameters of Fig. 15, for h = 1.091 and for h = 1.1. Although h = 1.091, which lies in the region 
of increasing magnetization (see Fig. 11), and is very close to the critical field ft.2c — 1-096, the 
oscillatory behaviour is still present in the correlation function. On the other hand, when h = 1.1, 
which Hes in a plateaux of the magnetization (see also Fig. 11), there is no oscillatory behaviour. 
It means that the period of the oscillation tends to infinity, and this result is still valid for any 
value of the field in the plateaux. This is consistent with the scaling form and the analytical 
continuation proposed for this correlation function [15], where the correlation length is associated 
to the oscillation period. 

Although the direct correlations (Sf^^S^^^J - {Sf^^) {S^+rJ and {Sf^Sf^.J - {SfJ (Sf^^J 
are not presented, they are easily obtained from the shown results and exhibit qualitatively the 
same behaviour as (•S'^^m'^j+r n) ~('5'/m) ('^j+r n)- This is a consequence of the fact that the 
critical fields are not dependent on the subcell. 

6 Conclusions 

We have considered the XY (s — 1/2) model on the alternating one-dimensional superlattice 
(closed chain), and an exact solution was obtained by using the Green function method. The exci- 
tation spectrum was determined, and explicit expressions were obtained at arbitrary temperature 
for the internal energy, the specific heat, the magnetization, the susceptibility, and the two-spin 
static correlation function in the field direction. 

The specific heat as a function of temperature, depending on the superlattice parameters, can 
present a single or a double peak, and we have shown that, at T = 0, dC/dT is different from zero 
provided there is a zero energy mode on the spectrum. 

In the T = limit, the behaviour of the system was studied as a function of the transverse field 
, and we have shown that the induced magnetization as a function of the field presents, alternately, 
regions of plateaus and of variable magnetization. Also in this temperature limit, the susceptibility 
in the field direction, x^^, presents singularities which are associated to phase transitions of second 
kind induced by the field. These critical points are consequence of the presence of zero energy 
modes with wave number or tt. In passing, it should be noted that this critical behaviuor, as 
expected, is suppresed at finite temperatures. 

These transitions have been treated within the real space renormaHzation group approach [16], 
and its critical exponents determined. The critical exponents can be also obtained directly from 
the exact expression, as shown elsewhere [13]. 

The two-spin static correlation function in the field direction, as a function of the separation 
between the spins, presents an oscillating behaviour in the regions where the magnetization is 
not constant, and the period of oscillations increases as the field approaches the critical value, 
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and diverges at h — he- This behaviour confirms that the static correlation function satisfies the 
analytical extension of the scaling form proposed for the homogeneous case [15]. 

Finally we would like to point out that the magnetization as a function of the field has qualita- 
tively the same behavior as those experimentally obtained for the NdCu2 in the low temperature 
limit, and these results have been obtained recently by Ellerby et al. [17] and Loewenhaupt et al. 
[18]. The agreement is more remarkable in the very low temperature limit, since in this limit the 
structure is analogous to a superlattice. Although this material is described by a Heisenberg type 
Hamiltonian, this result suggests that the lattice structure is a predominant factor in defining the 
magnetic properties of the material. 
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9 Figure Captions 

Fig. 1- Unit cell of tlie alternating superlattice. 

Fig. 2- Excitation spectrum of the homogeneous lattice for ua ~ '2,nB = 'S, J = Ja =, Jb = 1, 
and Ha = hs = h = 2. 

Fig. 3- Excitation spectrum for ha = 2, ub = 3, J = 1, Ja = 2, Jb = 3, and Ha = h-B = 0. 
Fig. 4- Excitation spectrum for ha = ns = 4, J = 1, J^ = 0.3, Jb = 3, and hA = h-B = 0. 
Fig. 5- Internal energy as a function of temperature (T* = fc^T) for n^ = n^ = 10, J = 1, 
Ja = 2, Jb = 3, and Ha — hB — 1-5. 

Fig. 6- Specific heat as a function of temperature (T* = fc^T) for nA = nB = 10, J = 1, Ja = 2, 
Jb — 3, and Ha — h-B ~ 1-5. 

Fig. 7- Specific heat as a function of temperature (T* = fesT) for ua = riB = 2, J = 1, Ja = 2, 
Jb ~ 3, and Ha — It-b — 1-3. 

Fig. 8- Specific heat as a function of temperature (T* = kBT) for ua = nB = 2, J = 1, Ja = 2, 
Jb = 3, and Ha ~ hB ~ 1.75. 

Fig. 9-The low temperature behaviour of the specific heat shown in Fig. 7 {T* = kBT). 
Fig. 10- The low temperature behaviour of the specific heat shown in Fig. 8 (T* = kBT). 
Fig. 11- (a) Average magnetization per subcell A and B (dashed line), and per unit cell (con- 
tinuous line) as a function of the field, and (b) the susceptibility in the field direction also as a 
function of the field, for ua = ub = 2, J ~ 1, Ja = 2, Jb = 3, and Ha = hB = h, at T = 0. The 
critical fields are hic = 0.691, /i2c — 1.096, h^c — 1.5960 and h^c = 1.809. 

Fig. 12- Average magnetization per unit cell as a function of the field, for ua — 2, hb = 3, J = 1, 
Ja = 2, Jb = 3, and hA ^ hB = h, aX T — 0. 

Fig. 13- Average magnetization per unit cell as a function of the field, for n^ = n^ = 4, J = 1, 
Ja — 2, Jb = 3, and Ha = 1.25/ib = h, at T = 0. 

Fig. 14- Average magnetization per unit cell as a function of the field, for ua = nB = 4, J = 1, 
Ja — 2,Jb ~ 3, and Ha ~ hB — h, at finite temperature (/3J — 50). 

Fig. 15- The direct static correlation function in the field direction, p^ (r) = {Sj^2^f+r.2) ~ 
{Sj^2 ) ; as a function of the distance between cells, for ua = tlb = 2, J = 1, Ja = 2, Jb = 3, and 
hA = hB^h, at T^ 0, for h = 0.75(a) and h = 0.693(b). 

Fig. 16- The direct static correlation function in the field direction, p"^ (r) = {Sj^2^f+r.2) ^ 
(3^2) ' as a function of the distance between cells, for n^ = n^ = 2, J = 1, J^ = 2, Jb — 3, and 
hA^hB^h,atT = 0, for h = 1.1(a) and h = 1.091(b). 
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